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; ABSTRACT 

Aims. A detailed study is presented, including estimates of the impact on elemental abundance analysis, of the non- 
(— I ' Local Thermodynamic Equilibrium (non-LTE) formation of the high-excitation neutral oxygen 777 nm triplet in model 

atmospheres representative of stars with spectral types F to K. 
I ■ Methods. We have applied the statistical equilibrium code MULTI to a number of plane-parallel MARCS atmospheric 

P ■ models covering late-type stars (4500 < T^s < 6500 K, 2 < log 3 < 5 [cgs], and -3.5 < [Fe/H] < 0). The atomic model 

employed includes, in particular, recent quantum-mechanical electron collision data. 
C/3 , Results. We confirm that the Ol triplet lines form under non-LTE conditions in late-type stars, suffering negative 

. abundance corrections with respect to LTE. At solar metallicity, the non-LTE effect, mainly attributed in previous 

studies to photon losses in the triplet itself, is also driven by an additional significant contribution from line opacity. At 
low metallicity, the very pronounced departures from LTE are due to overpopulation of the lower level (3s ^S°) of the 
transition. Large line opacity stems from triplet-quintet intersystem electron collisions, a form of coupling previously 
not considered or seriously underestimated. The non-LTE efi'ects become generally severe for models (both giants and 
dwarfs) with higher Tcff. Interestingly, in metal-poor turn-off' stars, the negative non-LTE abundance corrections tend 
to rapidly become more severe towards lower metallicity. When neglecting H collisions, they amount to as much as 
I A log eo| ~ 0.9 dex and ~ 1.2 dex, respectively at [Fe/H]= —3 and [Fe/H]= —3.5. Even when such collisions are 
included, the LTE abundance remains a serious overestimate, correspondingly by |Alog eol ~ 0.5 dex and ^ 0.9 dex at 
^SJ ' such low metallicities. Although the poorly known inelastic hydrogen collisions thus remain an important uncertainty, 

' the large metallicity-dependent non-LTE effects seem to point to a resulting "low" (compared to LTE) [O/Fe] in metal- 

^\ ^ poor halo stars. 

, Conclusions. 

Key words, line: formation - stars: abundances - stars: late-type - Galaxy: evolution 

C3 1. Introduction Type II supernovae (SNe II, Arnett 1978; Clayton 2003). 

Due to the longer lifetimes of their precursors, the delayed 

After the primordially ubiquitous H and He, oxygen is the release of additional iron from type la SNe starts to affect 

third most abundant chemical element in the Universe, and [q /FeQ only at a later stage. The oxygen abundance in un- 

the first among those produced via stellar nucleosynthe- evolved very metal-poor stars should therefore reflect the 

sis. It is a catalyst in the CNO process, which is the main enrichment by the very first generations of SNe. 

conversion channel of H into He in the interiors of stars The oxygen abundance at low metallicity has been the 

more massive than the Sun. The oxygen content of stars fo^us of a plethora of studies. Still, a long-standing debate 

and Its abundance ratio with other chemical elements de- grists on the exact [O/Fe] trend in metabpoor stars (see 

rived from observations of metal-poor halo stars are crucial g g_ discussion in Asplund 2005 and Melendez et al. 2006 

(e.g. Matteucci & Frangois 1992) to constraining Galactic ^nd references therein). Here, we wiU only briefly mention 

Chemical Evolution (GCE) models. Oxygen, like the other g^^^g pgj-tinent works related to the subject, 

a-elements, has been known for a long time (e.g. Conti et Neutral oxygen lines used in abundance studies of solar- 

al. 1967) to show an overabundance with respect to iron ^ypg ^^^^^ j^^j^^g ^-^^ excitation line at 615.8 nm, the 

m stars of low metallicity This is expected since its mam ^g^j^ forbidden [Oi] lines at 630.0 nm and 636.4 nm, and 

production site is m massive stars which end their lives as infrared triplet at 777 nm. The latter is the only strong 

oxygen atomic feature observable in spectra of late-type 
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stars. While it is known to suffer from large non-LTE ef- 
fects in early-type stars (Johnson, Milkey & Ramsey 1974; 
Baschek, Scholz & Sedlmayr 1977), there is little consen- 
sus on the estimated size of the abundance corrections for 
cooler stars of solar-type. For example, the formation of 
the IR triplet lines in the Sun has been the object of a 
number of investigations. With a semiempirical approach, 
Altrock (1968) found that the non-LTE source function 
drops below the Planck function in the solar photosphere, 
a result later confirmed by the multilevel non-LTE study by 
Sedlmayr (1974). More recently, the non-LTE abundance 
corrections in the Sun have been estimated to be between 

0.1 and -0.25 dex (Kiselman 1991, 1993; Gratton et al. 

1999; Takcda 2003; Takcda & Honda 2005; Asphmd ct al. 
2004), mostly depending on whether efficient collisions with 
neutral H atoms are adopted or not. These authors explain 
the non LTE line strengthening as due to photon losses in 
the lines, while the population (and thus the line opacity) 
of the lower level in the 777 nm transition is thought to 
stay close to LTE. 

Non-LTE effects on oxygen lines in solar-type stars 
have been reviewed by Kiselman (2001). Statistical equi- 
librium calculations using a 16-lcvel oxygen model atom 
(Kiselman 1991) have predicted an enhancement towards 
low metallicity of the (negative) non-LTE abundance cor- 
rections on the permitted 777 nm O l lines, suggesting that 
an LTE approach will seriously overestimate the derived 
oxygen abundance. The [0/Fe] ratio derived in LTE us- 
ing this abundance indicator would thus be increasingly re- 
duced in metal-poor stars once non-LTE effects are consid- 
ered. Other calculations have instead found almost [Fe/H]- 
indepcndent non-LTE corrections (Takeda 2003). 

Very high overabundances of oxygen compared to iron 
have occasionally been derived at low metallicity. Analysing 
the Oi 777 nm lines in metal-poor dwarfs, Abia & Rebolo 
(1989) foimd a steep monotonic increase in the oxygen-to- 
iron ratio for lower metallicities, reaching [0/Fe] > 1 below 
[Fe/H] ~ —2). Those authors investigated non-LTE effects 
for the triplet lines but claimed that they should be less 
severe than —0.2 dex, therefore still leaving a residual in- 
creasing trend. Their high abundance estimates were crit- 
icized in King (1993) as mainly due to too low effective 
temperature estimates. Other authors (Gratton & Ortolani 
1986; Barbuy 1988; Sneden et al. 1991; Spite & Spite 1991; 
Kraft et al. 1992: Carretta, Gratton & Sneden 2000: Cayrel 
et al. 2004; Spite et al. 2005) found [0/Fe]~ 0.4 - 0.6 in 
late- type stars with [Fe/H] < — 1, in particular using the 
[O l] lines in metal-poor giants. This nearly plateau-like be- 
haviour (with oxygen content varying much more slowly 
than iron at low [Fe/H], or even remaining constant) could 
be interpreted as the result of the contribution by type II 
SNe alone to the early Galactic nucleosynthesis (see e.g. 
Matteucci & Greggio 1986). While some authors have re- 
ported a fair agreement between the near-IR 777 nm oxygen 
lines and the forbidden line (Mishenina et al. 2000; Nissen et 
al. 2002), others (Takeda 2003; Fulbright & Johnson 2003) 
derived higher (up to several tenths of a dex) abundance 
values from the triplet. In particular even when applying 
to the observations the non LTE corrections derived in his 
own investigation Takeda still derived apparently too high 
(by several tenths of a dex) oxygen abundance from the IR 
triplet lines in metal-poor halo stars, a possible signal of 
imderestimated non LTE effects. At [Fc/H] < —3.5, even 
more extreme examples of the oxygen conflict have been 



found, with differences of up to 1.55 dex between [O i]-based 
and triplet abundances, and exceptionally large [0/Fe]> 2 
oxygen enhancements (e.g. Depagne et al. 2002; Israelian 
et al. 2004). 

Unfortunately, the neutral atomic oxygen lines, in par- 
ticular the LTE-obeying [O l] line, tend to become vanish- 
ingly weak in metal-poor solar-type stars. Thus, molecidar 
bands of OH in the UV (electronic transitions) and near IR 
(vibration-rotation lines) have also been used. The former 
are relatively strong and thus detectable in FGK stars down 
to low [Fe/H], the latter instead require high S/N and are 
only seen in spectra of cool stars (rcff< 5000 K). Results 
from both are however affected by likely important 3D ef- 
fects, in particular in the case of the UV lines (Asplund & 
Garcia Perez 2001), which may also suffer from the prob- 
lem of missing opacity (Balachandran & Bell 1998; Asplund 
2004) and possible issues related to continuum location, 
blends in the spectra, atmospheric extinction/cutoff and 
low CCD sensitivity. Various other uncertainties - e.g. on 
oscillator strengths for both molecular bands - and incom- 
plete knowledge of non-LTE effects on molecular line for- 
mation may also be important. Generally higher values 
of oxygen abundance (compared to those from the weak 
IR molecular features) are derived from the OH lines in 
the UV (Balachandran, Carr & Carney 2001; Melendez & 
Barbuy 2002; Barbuy et al. 2003), but the above effects 
may indeed already for the most part explain the large dis- 
crepancy. Other authors (Israelian, Garcia Lopez & Rebolo 
1998; Bocsgaard et al. 1999) have instead found reasonable 
agrcxanent between the O abundances derived for a sam- 
ple of mainly turn-off stars from the two indicators. In any 
case, it seems that consistent oxygen abundances using all 
the different diagnostics, namely permitted, as well as for- 
bidden and molecular lines, can hardly be obtained at low 
metallicity. Departures from homogeneity and from LTE 
have been at various times invoked for solving the inconsis- 
tency. A combination of other factors, including determina- 
tion of stellar parameters/temperature scale, quality of ob- 
servations and of reduced data, adopted solar abundances 
and so on may additionally be of concern and it is hard 
to disentagle the different effects. To complicate the issue 
further, one needs to consider that some authors have de- 
rived oxygen abundances at low metallicity from unevolved 
stars, others from giants, and it could be hypothesized that 
the lower oxygen values usually found for the latter are 
due to a depletion process in their atmospheres. However, 
Bessell, Sutherland & Ruan (1991) discussed evidence that 
oxygen abundances from OH were similar in dwarfs and gi- 
ants. More recently, Garcia Perez et al. (2006) confirmed a 
reasonable agreement between OH and [O i] abundances in 
mctal-poor subgiants, finding that O I triplet abundances 
give instead higher values even when corrected for non-LTE 
effects. It is therefore important to check if such residual 
discrepancy is linked to problems in the previous modelling 
of the formation of the permitted lines. 

Here, we thus focus on investigating departures from 

LTE in ID. Crucially, we explore the very low metallicity 
range, to understand if the non-LTE abundance corrections 
depend strongly on metallicity. In the next section, we de- 
scribe our non-LTE calculations, and give results in Section 
3. We then embark on a comparison with other non-LTE 
studies of the Oi777nm absorption lines (Sect. 4). Finally, 
the last two sections are dedicated to a discussion of the 



D. Fabbian et al.: Ol 777 nm triplet non-LTE line formation 



3 



2.2. Atomic model 




Fig. 1. Grotrian diagram of the atomic model we employed, 
showing the 53 energy levels of neutral oxygen and the sin- 
gle first ionization level. Lines connecting the various levels 
in the figure represent the 208 radiative bound-bound tran- 
sitions included in the model. The 777 nm transitions are 
marked with thicker lines 



impact on the [0/Fe] ratio at low metallicity, and to some 
concluding remarks. 



2. Non-LTE calculations 

2.1. Method 

The method employed is similar to that described in 
Fabbian et al. (2006). In order to derive the non-LTE pop- 
ulations of the O atomic levels and the strength of rel- 
evant spectral lines, it is necessary to solve the coupled 
rate and radiative transfer equations simultaneously. The 
code MULTI (Carlsson 1986) in its version 2.3 was used to 
perform the statistical equilibrium calculations for a num- 
ber of MARCS model atmospheres (Gustafsson et al. 1975; 
Asplund et al. 1997), under the assumption that oxygen is 
a trace element and its departures from LTE will not affect 
the atmospheric structure. The atmospheric models tested 
have a microturbulence of ^ = 1. A 54-level oxygen model 
atom (described below) was employed. We assumed com- 
plete redistribution in a Voigt profile for the formation of 
the lines. 



We have constructed a model atom containing 54 en- 
ergy levels and a total of 258 radiative transitions (208 b-b 
and 50 b-f). The Grotrian diagram of our atomic model is 
shown in Fig.[T] The necessary data for atomic level ener- 
gies and corresponding statistical weights were taken from 
the NIST Atomic Spectra DatabaseQ version 3. Neutral lev- 
els are included up to an excitation potential of 13.34 eV 
(i.e. ~ 0.28 eV below the continuum). The model is com- 
plete to a principal quantum number of n= 7. Fine-splitting 
of energy levels was taken into account where appropriate, 
in particular for the ground state and for the 3p level 
(upper level of the 777 nm Oi triplet). For the transitions 
of importance in solar work which were used in Asplund 
et al. (2004), f- values are also from NIST while radiative 
and Stark parameters come from VALE0 (Vienna Atomic 
Line Database, Piskunov et al. 1995 and successive up- 
dates). Data for bound- free photoionization are from the 
Opacity ProjeclQ (Seaton 1987 and later updates). For col- 
lisional broadening due to H atoms, we adopt the classical 
Unsold approximation. Our tests reveal that choosing an 
enhancement factor of two for the related damping con- 
stant does not affect our results significantly. The contribu- 
tion to background opacity in the UV from lines of other 
elements was accounted for by compiling NIST and VALD 
values. In general, the data in the model atom^ of Carlsson 
& Judge (1993, hereafter CJ93) and Kiselman (1993, here- 
after K93) were used (see those papers for description of 
the data sources) unless otherwise stated. 

We have also included the most recent data for electron- 
impact excitation, available up to level 2p'^ 4f ^F from 
quantum- mechanical calculations (Barklem 2007a) , a novel 
feature of our model atom which has proven one of the 
essential improvements with respect to other studies. In 
particular, this has allowed us to include accurate rates for 
transitions between radiatively forbidden transitions due to 
electron collisions. 

Since the estimated cross-sections were derived using 
LS-coupling, they are null only for singlet-quintet transi- 
tions (total electronic spin conservation rule) and for coUi- 
sionally forbidden S'^ - S° transitions. In broad terms, the 
collisional cross-sections compare fairly well with available 
empirical data (within the experimental error bars, see Fig. 
2 in Barklem 2007a). Although admittedly still subject to 
uncertainties due to e.g. pseudo-resonances at higher ener- 
gies, they should be more realistic than previous estimates. 
They are in general comparatively larger than H collisions 
by a few orders of magnitude, except for fine-splitting lev- 
el^. We stress here that the cross-sections for electron- 
impact collisions in the radiatively forbidden transition 3s 
"^3° - 3s ^S° are among the largest in the calculations of 
Barklem (2007a). In particular (see Fig. 1 in that work) 
they are larger than the other values for transitions arising 
from the lowest levels (and between the first seven energy 
levels in our atomic model). 

^ http:/ /physics. nist.gov/PhysRefData/ASD/index. html 
^ Available at "http:/ /ams. astro. univie.ac.at/vald/| 
* http:/ /vizier. u-strasbg.fr/topbase/topbase. html 
^ Respectively, a 45- and a 23-level atomic model 
^ We account for fine splitting also for electron-impact data, 
by using Barklem's calculations (which originally refer to 
grouped levels) according to conservation of total rates 
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Due to the deficiency of metals, and thus the decreased 
number of free electrons in their atmosphere, rates due 
to collisions with neutral H atoms may be particularly 
important in low-metallicity stars. The approximation we 
adopted for excitation and ionization via inelastic collisions 
with Hi is based on recipes by Steenbock & Holweger 
(1984), which further generalize results in Drawin (1968, 
1969), where Thomson's classical theory for electron-atom 
encounters was applied for collisions between identical par- 
ticles. The approach is admittedly less than ideal, but un- 
avoidable due to the lack of relevant experimental data and 
theoretical calculations (see e.g. the discussion in Asplund 
2005). In fact, quantum mechanical calculations of rate 
coefficients exist for only very few electronic transitions 
of oxygen atoms induced by collisions with H (Krems, 
Jamieson & Dalgarno 2006; Abrahamsson et al. 2007), and 
they are unfortunately often of limited applicability for stel- 
lar work due to the low range of temperatures explored. 

Thus, in order to estimate the size of the main uncer- 
tainties in the calculations, we have carried out, similarly to 
what was done in Fabbian et al. (2006), different non-LTE 
calculations with a scaling factor Sh for transitions between 
all levels, including for radiatively forbidden transitions (by 
adopting a minimum f- value of 10^^ in that case). 

3. Results 

3.1. Non-LTE mechanisms 

The behaviour of the departure coefficients /3i (defined 
as the ratio of the populations in non-LTE and LTE, f3i = 
rii/nY^^) with atmospheric depth for various energy levels 
in the atomic model is shown in Fig. [5] for the Sun and a 
metal-poor turn-off model. 

Non-LTE effects for the O i 777 nm lines at solar metal- 
licity are mainly caused by source function dilution, in the 
atmospheric layers where the lines form, with respect to 
the Planck function that represents the LTE expectation. 
This is due to radiative losses in the triplet itself causing 
line strengthening compared to LTE. As seen in Fig. [21 the 
upper level (3p ^P) of the transition gets underpopulated, 
while the lower level (3s ^S°) gets overpopulated. If stim- 
ulated emission is neglected, this driving effect can be well 
described as Si/B = Pu/Pi- The two-level approximation 
will give a reasonable description of the non-LTE line for- 
mation in the Sun, implying that it is the processes in the 
transition itself that matter most. 

At low metallicity instead, particularly at the highest 
temperatures explored, that approximation breaks down. 
The upper level (3s '^S°) of the 130 nm Oi resonance lines 
is in effect metastable, since radiative de-excitation will be 
immediately followed by an excitation (due to large opacity 
in those lines). Thus, very large overpopulation compared 
to LTE develops in the level, because of the larger net rate 
for it than in LTE. This is caused by radiative transitions 
having this as their lower level (i.e. by photon losses in 
higher excitation levels), combined with the low rate coef- 
ficients out of the level. 

The 3s ^S° overpopulation tends to propagate to the 
3s ^S° level via efficient intersystem collisional coupling, 
and line strengthening of the 777 nm lines results, due to 
much increased line opacity. Given the large energy differ- 
ence with the ground state, collisions to it from the levels 
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Fig. 2. Departure coefficients in the level population of the 
lowest energy levels (up to 3p ^P) in our oxygen model atom 
(for the case without H collisions) are shown respectively 
for the Sun (upper panel) and a metal-poor turn-off at- 
mospheric model with [Fe/II]= —3 (lower panel). Note the 
very different scales for the y-axes in the two panels. The 
thicker lines indicate the departure coefficients for lower 
and upper levels (solid and dashed curve respectively) in 
the 777 nm triplet. Note that only one curve is plotted for 
the latter levels, because the small differences in the depar- 
ture coefficients among those fine structure levels does not 
play an important role here 



of interest are inefficient in maintaining LTE. The domi- 
nant contribution to the non-LTE effect at low metallicity 
thus comes from line opacity, because the lower level of the 
IR triplet is strongly overpopulated with respect to LTE 
(Fig. (2). The non-LTE abundance corrections then become 
^ —0.5 dex at [Fe/II]= —3 for a typical metal-poor turn-off 
halo star, and up to —0.9 dex if neglecting H collisions. 

As seen in Fig.[3l in non-LTE the line formation at 
very low metallicity is shifted to higher atmospheric lay- 
ers. Since the population of the lower level in the oxygen 
777 nm triplet transition rises very steeply with respect to 
the LTE expectation when moving outward in the atmo- 
sphere, the P departure coefficient will become large. Fig. [3] 
also shows that the effect (on line-center flux) of photon 
losses in the line, causing the source function to drop below 
the Planck function, is small at low [Fe/H]. The opacity 
effect dominates, with optical depth unity moving up in 
the atmosphere because of the overpopulation of the lower 
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level of the transition. The large line opacity at the shal- 
lower depths where line formation occurs in non-LTE will 
thus make the absorption stronger. 

Even with decreased free electrons at low metallicity, 
sufficient intersystem coUisional coupling is maintained be- 
tween levels of similar excitation in the triplet and quintet 
systems. The 3s '^S° and the 3s levels are in fact very 
close energetically and, despite the corresponding transition 
being radiatively forbidden, they are strongly coupled via 
collisions, which explains why their departure coefficients 
behave similarly, as seen in Fig. [21 Due to decreased num- 
ber of electrons and to small rates for coUisional transitions 
with larger energy separation, electron collisions are more 
and more ineffective in driving the line formation closer to 
LTE via depopulating these two levels via other channels, 
e.g. to the ground state. The less efficient thermalization 
via impacts is thus unable to balance the tendency of the 
triplet levels to overpopulate. In addition, the lower level 
of the 777 nm transition is effectively metastable due to 
negligible net radiative rates to the ground state. 

Among the IR Oi triplet lines, the strongest (bluest) 
suffers generally more severe non-LTE effects. This is due 
to it having relatively more important photon losses, and 
given that it will form further out where departures from 
LTE are largest. The difference in the estimated abundance 
correction among the 777 nm lines usually remains within 
0.05 dex, but grows towards solar metallicity. The abun- 
dance corrections differ by as much as ^ —0.15 dex at 
Toff= 6500 K, log g= 2, [Fe/H]= 0, where very large non- 
LTE effects are found. Thus, adopting the same correction 
for all three lines is not appropriate, especially in such case, 
and detailed non-LTE calculations are necessary. This is 
due to substantially different line formation depths, giving 
larger non-LTE effects for the line that is formed further 
out, i.e. the strongest (bluest) one, because of outward in- 
crease of the departure coefficients. 

Across all of the parameter space explored, the upper 
level (3s ^S°) in the Ol UV resonance lines tends to get 
overpopulated. Large radiation excess is maintained in the 
resonance lines, becoming particularly strong at low metal- 
licity, for high gravity models. 

In the 777 nm lines instead, the ratio of J over B tends 
to be < 1 due to photon losses. This is relatively more 
important in stars like the Sun, and causes the upper levels 
in the 777 nm triplet to depopulate. 

Photoionization was found not to be important across 
the parameter space explored in this study. 



Table 1. Results (for the strongest triplet line at 777.2 nm) 
obtained using our 54-level atomic model. The non-LTE 
and LTE equivalent widths are given to show the sensi- 
tivity to various mechanisms. The changes with respect to 
the standard case (Teff= 6500 K, log g= 4, [Fe/H]= -3, 
^ =1, log eo = 6.5, and including the coUisional data of 
Barklem 2007a) are: decrease in metallicity by half dex (to 
[Fe/II]= —3.5); and decrease in oxygen abundance by half 
dex (to log eo = 6.00). For the case including H colUsions, 
additional tests are, following the order listed: removing 
the 130 nm O I resonance lines; removing also the intersys- 
tem 136 nm UV O i lines; removing all of the new electron 
collisions by Barklem except in the radiatively forbidden 
transition between levels 3s '^S° and 3s ^S°; including all 
the new electron collisions by Barklem except the previ- 
ous one; no 3s and 3s coupUng (i.e. removing both 
electron and H collisions between those levels); use of old 
electron collision data from 23-level model atom of K93 in- 
stead of Barklem's data; same as previous but with new 
electron collision data between levels 3s '^S° and 3s ^S°; old 
electron data but no H collision coupling between the latter 
levels; neglect of UV background line opacities; exclusion of 
the opacity contribution by the blending H Ly-/3 line; and 
decrease by 1 dex in Si abundance 
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[mA] 
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[mA] 
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standard case 


14.2 


2.9 


[Fe/H]= -3.5 


10.9 


1.0 


log eo = 6.0 


6.3 


1.0 



Sh = 1 



standard case 


6.7 


2.9 


[Fe/H]= -3.5 


4.2 


1.0 


log eo = 6.0 


2.9 


1.0 


no 130 nm O I lines 


3.9 


2.9 


no 130 nm, no 136 nm O I lines 


3.6 


2.9 


3s - 3s Barklem e~ 


7.2 


2.9 


all Barklem e~ but 3s - 3s 


6.5 


2.9 


no 3s - 3s coupling 


5.2 


2.9 


old e~ collisions 


6.2 


2.9 


old e~ collisions, 3s - 3s Barklem e~ 


6.4 


2.9 


old e~ collisions, no O+H 3s - 3s 


5.0 


2.9 


no UV background line opacity 


7.8 


2.9 


no H Ly-/3 opacity 


6.8 


2.9 


Si*0.1 


9.4 


2.9 



3.1.1. Sensitivity of line formation: main atomic transitions 



In Table [T] we provide the equivalent widths resulting 
from some of the various tests carried out on our 54-level 
atomic model in order to check the sensitivity of the 777 nm 
triplet line formation to various transitions and processes in 
the atomic model. The standard case shown (Toff = 6500 K, 
log g= 4, [Fe/H]= -3, ^ =1, log eo = 6.5) includes the 
new electron-impact data by Barklem (2007a). It results 
in A log eo ~ —0.9 and ^ —0.5 dex without and with 
the inclusion of H collisions respectively. The correspond- 
ing non-LTE abundance corrections at [Fe/II]= —3.5 in- 
crease to A log eo ~ —1.2 and ~ —0.9 dex respectively. 



The modifications were performed on this standard case, 
as listed in Table [TJ From the results there shown, it is 
clear that the most important transitions are the O i res- 
onance lines (together with the corresponding UV back- 
ground opacity), and the intersystem coUisional coupling 
between the 3s ^S° and 3s ^S° states. The model atom em- 
ployed in the standard case is also the one we then used to 
derive our non-LTE corrections across the whole parameter 
space (see Section lX^ . 

We also carried out "multi-MULTI" tests, with results 
for the case without H collisions shown in Fig. [4] (only avail- 
able online). Through such tests, the relative effect on the 
O I 777 nm triplet line strength of multiplying each of the 
different rates for radiative and electron coUisional transi- 
tions individually by a factor of two is revealed. One can 
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Fig. 3. The various panels refer to the oxygen feature at 777.2 nm (the bluest and strongest of the Oi triplet lines), for 
[Fe/H]= —3.00 (top panels) and —3.50 (bottom panels). The atmospheric models have Tcs— 6500 K, log g— 4.00 [cgs] 
and solar-scaled, alpha-enhanced ([0/Fe]= -1-0.4) oxygen abundance. The case with no hydrogen collisions is shown here. 
Top: the resulting LTE (dashed curve) and non-LTE (solid curve) flux profiles, normalized to continuum flux, for the 
spectral feature considered are shown in the left panel, while normalized values of B (dashed), S (solid) and Si (dash-dot) 
are shown in the right panel. The vertical scale in the two cases is therefore the same. The left panel illustrates how at 
low metallicity the line becomes much stronger in non-LTE compared to LTE. In the right panel, vertical lines show the 
continuum (dotted), LTE line center (dashed) and non-LTE line center (solid) Eddington-Barbier values of for which 
the continuum flux corresponds to the Planck function at that depth (i.e. F™"*=7rSiy(T™"*=E.B.)). The corresponding 
starred symbols on the various curves therefore represent the Eddington-Barbier formation height and flux, respectively 
for continuum, LTE line center and non-LTE line center. Note that the Eddington-Barbier value is r^°"* = 2/3 exactly 
when the source function varies linearly with optical depth, while it is slightly smaller in this case. The right panel shows 
that at low metallicity the source function effect is small. Bottom: as above, but for [Fe/H]= —3.50 



thus grasp the impact of different atomic transitions on the 
formation of the lines of interest. As seen in Fig. 01 and 
as found by previous authors (e.g. Kiselman 1991, 1993), 
the formation of the O i triplet in the Sun is mostly influ- 
enced by photon losses in those lines and can be described 
reasonably well by a two-level approximation. We stress 
however that our calculations also show an additional con- 



tribution from line opacity, due to the lower level of the 
triplet transition getting overpopulated. This increases the 
resulting non-LTE effects for the Sun to a larger amount 
than usually adopted, up to —0.3 dex when neglecting H 
collisions. 

At these solar or moderately low metallicities, the line 
opacity effect in the O i IR triplet starts to dominate (over 



D. Fabbian et al.: Ol 777 nm triplet non-LTE line formation 



7 



Fig. 4. (only available online) Results from multi-MULTI runs for a solar model, and for a star with T^s— 6500 K, 
log g= 4.00 and [Fe/H]= —3, when neglecting H collisions. 



that caused by the source function drop) for hot giants. It 
becomes much more predominant at very low metallicity, 
where we find an increase in the size of the abundance cor- 
rections with increasing Tcs and decreasing [Fe/H], related, 
respectively, to increased level overpopulation via collisions, 
and to decreasing continuum opacity. An increase in gravity 
at metallicities below [Fe/H]'--^ —2 will also tend to cause 
larger non-LTE corrections. This effect is related to the 
increased sensitivity of the line formation process to inter- 
system collisional coupling, in that regime. 

So, while in the Sun both source function drop and line 
opacity are important, with the first being predominant, 
the metal-poor turn-off model is mostly affected by the 
latter. The effect of the line source function dilution does 
not change significantly with [Fe/H] at very low metallic- 
ity, its contribution to line strengthening remaining < 10%. 
For parameters typical of metal-poor turn-off stars, we find 
a trend of increasing non-LTE corrections with decreas- 
ing [Fe/H]. This is caused by the behaviour of continuous 
opacity at the UV wavelengths of interest (i.e., around the 
Oi 130 nm features), which decreases substantially (up to 
~ 15%) across the line formation layers when the metal con- 
tent goes from [Fe/H]= —3 to —3.5. Even though the con- 
tribution of Rayleigh scattering to total continuous opacity 
is dominant for layers around and above those where the 
777 nm features are formed, it does not vary when metallic- 
ity decreases from [Fe/H]= —3 to —3.5. Instead, it is the de- 
creased contribution by Si photoionization to total continu- 
ous opacity when [Fe/H] decreases (see Fig. [5]) that deprives 
the photons in the O i resonance lines of an alternative 
route. Their destruction probability via the photoionization 
channel therefore decreases. As a consequence, overpopula- 
tion (driven by photon losses in higher-excitation lines) of 
the upper level in the O i resonance lines tends to increase. 
At the same time, the lower level of the 777 nm transition 
will too correspondingly get more overpopulated due to ef- 
ficient intersystem collisional coupling. At low metallicity, 
the non-LTE effect on the O i 777 nm triplet feeds on the 
flow via collisions from the 3s state. The level occu- 
pation numbers in the lower level of the IR oxygen triplet 
thus stay closely coupled to the large overpopulation in the 
upper level of the UV 130 nm resonance lines, a situation 
which LTE cannot predict. 

At [Fe/H]'--^ —3, the triplet lines become stronger when 
increasing the radiative rates for the resonance lines, be- 
cause that causes larger overpopulation in level 3s '^S°, and 
in the lower level of the 777 nm triplet as a consequence 
of collisional coupling. Note that in the figure, we evaluate 
the net flow (indicated by the direction of the arrows for 
each transition) where r = 1 at line center for the 777 nm 
triplet. This generally gives a good indication of the driv- 
ing mechanisms. However, at low metallicity, one can note 
that the radiative rate in the triplet is upward. A check re- 
veals that in this case the rates actually change sign across 
the atmosphere. While at r = 1.0 at line center there is a 
small upward flow, the net rate turns to being downward in 
deep layers. This last effect dominates, explaining the over- 
all line strengthening, as marked in the corresponding fig- 
ure, when increasing the radiative rate in the O i IR triplet. 



Excluding the transition itself, the most significant impact 
on the triplet lines at low metallicity is from pumping in 
the Oi 130 nm resonance lines. 

Concerning electron collisions, the main effects at low 
metallicity are strengthening of the line through relatively 
efficient triplet-quintet intersystem coupling via electron 
collision^ between the similarly highly-excited levels 3s 
and 3s ^S° - which helps to propagate the overpopula- 
tion from the former state to the IR oxygen triplet - and 
line weakening via collisions to the ground state. At low 
metallicity, both lower (3s ^S°) and upper (3p ^P) levels 
of the triplet lines get very overpopulated, but by different 
amounts. 

3.1.2. Importance of Oi UV transitions and background line 
opacities 

According to our calculations, even in metal-poor stars, 
at least for frequencies near their quasi-saturated core, 
the extremely optically thick 130 nm Oi resonance lines 
(2p'* ^P2,i,o-3s 3S°) form in the very topmost layers of 
the photosphere as strong absorption features. Despite this, 
their thermalization depth lies much deeper in than where 
monochromatic optical depth in the lines is unity, due to 
the very large amount of scattering. In particular their 
J/B>> 1 where the 777 nm lines are forming. The prop- 
erties of scattered photons are decoupled from local condi- 
tions, which gives a typical non-LTE situation. 

The fact that the O i resonance lines are so optically 
thick even at low [Fe/H] means that at the formation lay- 
ers of interest, collisional de-excitation from their upper 
(metastable) level to the lower level of the 777 nm transi- 
tion will be favoured with respect to de-excitation to the 
ground. This situation causes large non-LTE effects at low 
metallicity. When removing the UV O i resonance lines, 
the non-LTE effects remain important (see Tabled]) due to 
photon losses from higher levels. 

Previous studies (e.g., K93; Gratton et al. 1999) found 
that UV radiation should not play a major role in the for- 
mation of the IR oxygen triplet in very-metal poor star. 
They were, however, lacking the detailed quantum mechan- 
ical calculations for intersystem electron collisions which we 
employ here and which cause increased non-LTE effects. 

Background absorption from UV lines of other elements 
mitigates (by providing an escape route other than col- 
lisional de-excitation and collisional transfer to the quin- 
tet system) but does not extinguish this effect. The ex- 
act amount of bound-free absorption from different met- 
als is however still debated (e.g. Balachandran & Bell 

^ As seen in Fig.|4l this channel gives instead a weakening ef- 
fect in the Sun. The reason is clear from the previously discussed 
Fig. [2] which shows that in the Sun, the departure coefficient is 
larger for 3s ^S° (lower level of 777 nm transition, in the quintet 
system) than for 3s '^S° (upper level of O I 130 nm UV reso- 
nance lines, in the triplet system). An increase of this collisional 
coupling channel in the Sun will therefore decrease the overpop- 
ulation and thus the opacity effect for the 777 nm transitions, 
making those lines form deeper in, where smaller non-LTE ef- 
fects are felt. 
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Fig. 5. Contributions to total continuous opacity around 
the 130 nm O i lines for a Tes= 6500, log g= 4, [Fe/H]= -3 
model atmosphere. The absorption and the scattering part 
of total opacity are shown by the two thick dashed curves, 
as marked. For absorption, the importance of its main com- 
ponents (Si, H^, and H I) is depicted. Total continuous 
opacity itself is shown by the thick solid curve. For compar- 
ison, we also show (small circles) the reduced contribution 
of Si at [Fe/H]=: —3.5 and consequent reduced total con- 
tinuous opacity (large circles). The rest of the sources of 
opacity in the figure do not change significantly with this 
decrease in metallicity 



1998; Asplund 2004), due to the possibility that contin- 
uum rates, e.g. for Fe I, may be larger than what usually 
adopted. As noted by AUende Prieto, Hubeny & Lambert 
(2003), predicting a realistic flux in the UV is more compli- 
cated than for visible or IR wavelengths. The knowledge of 
accurate bound-bound and photoionization cross-sections 
and thresholds for metals becomes important, with a large 
quantity of accurate atomic data required for modelling. 
Below ~ 1700 A, Si I bound-free absorption plays an im- 
portant role. Rayleigh scattering by hydrogen atoms is also 
important in the ultraviolet for cool stars, in particular at 
low metallicity. In addition to bound-free metal edges, a 
sharp increase in opacity is expected at short UV wave- 
lengths also due to the haze caused by background line 
absorption. Neglecting these effects will seriously overesti- 
mate the UV radiation field within the upper atmospheric 
layers of late-type stars. 



Fig. [5] shows the main contributions to total continu- 
ous opacity for UV wavelengths around the O i resonance 
lines and a MARCS model having Tcff^ 6500 K, log g= 4, 
[Fe/H]= —3 and log 0= 6.5. As seen, scattering is roughly 
constant across the atmospheric model. Its contribution 
to total opacity dominates down to the atmospheric lay- 
ers (log Tsoo ~ —0.5) where the IR oxygen triplet lines 
form. The contribution of absorption, even though smaller, 
is still significant, at such layers predominantly thanks to 
Si I bound- free edges. It is this component which proves 
important to understand the trends of the non-LTE cor- 
rections at low metallicity, since it is the only one which 
varies significantly with [Fe/H]. At [Fe/H]= —3.5, its con- 
tribution to absorption for the atmospheric layers of inter- 
est reduces to ~ 35% of the value at [Fe/H]= —3. Thus, 
total continuous opacity becomes significantly smaller. No 
other metal component plays a significant role. Deeper in, 
absorption takes over, due to and at first, and then 
to the rapid increase of neutral hydrogen photoionization. 

The larger continuum absorption at [Fe/H]= —3 com- 
pared to [Fe/II]= —3.5 reduces the overpopulation in the 
upper state of the UV Oi resonance transitions because 
it provides an alternative route to photons in those lines 
(absorption by the continuum instead of by the line itself). 

The non-LTE corrections are found to be very sensitive 
to the UV radiation field. Thus, background opacity caused 
by lines of other elements has been accounted for in the case 
of important transitions, in particular the contribution of H 
Ly-/3 and of features around the O i 130 nm resonance lines. 
The latter inclusion has proven particularly important be- 
cause it significantly reduces the otherwise even more ex- 
treme non-LTE effects (see Table [ij. In total, 151 lines 
from VALD were used, around the 130, 103 and 136 nm 
transitionfQ. For these background features, the line source 
function was assumed to be equal to the Planck function. 
This approximation is sufficient also to describe the case 
of pumping by hydrogen Ly-/3 radiation. This was tested 
by solving the radiation transfer problem for the hydrogen 
atom and taking the J field as input for the non-LTE cal- 
culations on oxygen. The results showed only very small 
differences with respect to our approximation, for the for- 
mation of the oxygen triplet lines of interest here. Note that 
H Ly-/? gives the single largest contribution to background 
line opacity, however it is the total contribution of the sev- 
eral blending metal lines falling around 130 nm which is 
most important, in particular due to Si 130.23 nm, which 
falls close to line center of the bluest O i resonance line. 
Sin 130.44 nm, which is located in the blue wing of the 
central resonance line but also affects the 103.22 nm Oi 
line by depressing the continuum, and S i 130.59 nm, which 
falls close to the center of the reddest resonance line. In 
addition, a number of lines of Fe, Ca, P and other elements 
give a smaller contribution to line haze. The lines around 
the Oi UV intersystem lines are instead not important. 

Our study reveals that the inclusion of the UV opac- 
ity contribution does not have a significant impact on the 
777 nm triplet lines in the Sun, where the non-LTE effect 
is in that case at least in good part controlled by processes 
in those same IR lines. 



* Respectively accounting for, in order of decreasing impor- 
tance, the effect of features around the resonance lines, of ab- 
sorption due to H Ly-/3, and of lines falling close to the inter- 
system oxygen lines 
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Fig. 6. Net rates contributing to the balance of levels 3s 
(left) and 3s (right). Lines with and without diamond 
symbols respectively indicate the contribution from colli- 
sions and radiation to a given rate. Only the rates with large 
influence on the level population are shown (their summed 
contribution is represented by crosses). The thick horizontal 
line at dn/dt= marks the statistical equilibrium expecta- 
tion that the total sum of all rates to and from a given level 
is equal to zero. The atmospheric model has T^b^ 6500, 
log g= 4 and log eo = 6.5. Plots are for [Fe/H]= —3 (up- 
per panels) and —3.5 (lower panels) respectively 



3.1.3. Effects of collisions with electrons and hydrogen 

Figure [H] shows the net rates at low metallicity for the 
high-excitation levels 3s '^S" and 3s ^8°. The fi gure clearly 
illustrates how the balance in the population of the lower 
level of the 777 nm oxygen triplet is mainly determined by 
incoming flow due to coUisional coupling - which transmits 
the overpopulation present in the 3s state - moderated 
by radiative excitation in the IR lines. 

We have included in our model atom a homogeneous 
set of electron-impact excitation data, thanks to the avail- 
ability of recent theoretical quantum mechanical calcula- 
tions (Barklem 2007a) . Table [T] shows the results of testing 
various modifications to the atomic model, including for 
electron collisions. The effect of the new electron-impact 
data is to give larger non-LTE corrections due to increased 
intersystem coupling, in particular, in the 3s - 3s 



transition. This is an efficient channel due to its large cross- 
section (it is a spin-ffip transition induced by an exchange 
interaction of the electrons) and its small 0.38 eV) en- 
ergy separation (thus low threshold). 

Inelastic collisions with neutral H atoms as modelled 
here can play a crucial role at low metallicity, where the 
scarcity of free electrons means that rates due to H colli- 
sions can become important in late-type stars due to the 
large H population density. Under the assumption that van 
Regemorter's and Drawin's formulae well reproduce the 
thermalization due to electron and H collisions respectively, 
one derives that the ratio of the rates due to the two pro- 
cesses is a function of the energy separation between the 
levels involved and of temperature ( Asplund 2005) . These 
collisions can thus become particularly efffcient for levels 
that are very close in energy. The approximate recipe avail- 
able is based on work by Drawin (1968, 1969) which general- 
izes the modified classical Thomson formula to atom-atom 
collisions. Unfortunately, as already observed by Steenbock 
and Holweger, Drawin's formula allows at best an order- 
of-magnitude estimate of the importance of H collisions. 
Significant differences have been found compared to more 
reliable estimates (e.g. experimental data by Fleck et al. 
1991; calculations for Na by Belyaev et al. 1999; for Li by 
Belyaev & Barklem 2003; and for H-hH by Barklem 2007b), 
differing by up to six orders of magnitude for some transi- 
tions. The uncertainty on the efficiency of H collisions has 
a significant impact on non-LTE studies. There is some ev- 
idence (Allende Prieto, Asplund & Fabiani Bendicho 2004) 
that for oxygen, relatively efficient collisions with H need 
to be included in statistical equilibrium calculations. We 
confirm the impact of their thermalizing effect, especially 
at low [Fe/H] (see Table [1]). They tend to mitigate the large 
overpopulation in the lower level of the 777 nm triplet and 
thus reduce the line opacity. The non-LTE abundance cor- 
rections at low metallicity become considerably smaller (by 
> 50%). 

In the absence of more reliable calculations, the only 
possibility to have some indication of the inffuence of H col- 
lisions has been to adopt the Drawin recipe. Even though 
likely not realistic, this may be useful in the case of oxy- 
gen, at least (inevitably) as a first approach, to investi- 
gate the impact of H collision efficiency. Given available 
evidence pointing towards the Drawin formula generally 
overestimating rates for individual transitions, the range 
of Sh = — 1 adopted here would seem reasonable for our 
tests, and may be expected to give an idea of the uncertain- 
ties involved with respect to stellar abundance estimates 
from the 777 nm Oi triplet. 

We also tested the effect of including the recent data 
by Krems, Jamieson & Dalgarno (2006). Even though their 
estimated rates (using formula 10 in their paper) are larger 
than our adopted values, they did not have a significant 
impact on our results since the transition they consider in 
the singlet and triplet systems (^D-'^P) involves levels that 
follow LTE anyway. 

3.2. Overall non-LTE effects across the parameter space 

The high-excitation triplet lines are particularly sensi- 
tive to temperature and become much stronger, both in 
LTE and non-LTE, for larger T^s values, because the cor- 
responding level population increases, i.e. there are more 
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Fig. 7. Non-LTE abundance corrections versus metallicity, 
for the strongest O i IR 777 nm triplet line. The lines con- 
nect models with given TcS and log g but varying metal- 
licity. Different symbols denote, respectively, results for 
models representative of a normal dwarf {Teg= 5780 K, 
log g= 4.44, indicated by open circles), a normal turn-off 
star (Teff= 6500 K, log g— 4, marked with filled diamonds) 
and a RR-Lyrae star (Toff = 6500 K, log g= 2, indicated by 
open squares). Dashed lines are for no H collisions, solid 
lines for a choice of Sh = 1- Solar-scaled oxygen abundance 
with a-enhancement below solar metallicity was adopted 
in the calculations, as indicated by the horizontal axis at 
the top of the figure 



atoms with electrons in the proper energy level to produce 
the spectral lines. Surface gravity (pressure) controls the 
amount of recombination via the Saha equation, but also 
the relative importance of coUisional processes with respect 
to radiative ones. The lines get stronger in lower gravity 
models. Metal content - i.e. smaller or larger number of 
absorbers - also controls the strength of the lines. 

An IDL routin43 has been prepared in order to interpo- 
late among the grid, therefore making it possible to calcu- 
late detailed abundance corrections for any given star with 
parameters in the range covered. 

The non-LTE line strength is larger than in LTE across 
the parameter space explored, i.e., negative abundance 
corrections arc found. Fig. [7| shows in what measure the 
strongest (bluest) Oi triplet line at 777.19 nm is affected 
by the resulting non-LTE corrections for some of the mod- 
els (typical dwarf, typical turn-off star, and RR Lyrae star). 
A qualitatively similar trend is seen, namely of decreasing 
abundance corrections for all models from solar metallicity 
to [Fe/H]= —2 and then (except when including H colli- 
sions in the RR Lyrae case) of increasingly severe non- 
LTE corrections towards very low metallicity, in particu- 
lar for turn-off stars, where the resulting |Alog eo| reaches 
~ —0.9 dex at [Fe/H]= —3, when neglecting H collisions. 
This reduces to ~ —0.5 dex if H collisions are included. 
For the atmospheric models considered, non-LTE correc- 

^ The routine and associated data are made available for gen- 
eral use. They are downloadable either as individual files, or as 
a single compressed file, from ftp:/ /nlte:cl olnlte@ftp. iac.es (at 
this same address the routine for carbon non-L'i'E corrections 
from Fabbian et al. (2006) is also available) 



tions are less severe around [Fe/H]= —2, but still remain 
significant, reaching ^ —0.4 dex for the RR Lyrae star. 
Interestingly, especially for the latter, H collisions do not 
have much impact (the difference staying within 0.1 dex 
for all models) unless [Fe/H]< —2. At higher metallicity, 
including H collisions in the RR Lyrae case actually makes 
the non-LTE abundance corrections ever so slightly more 
severf°l. For the Sun, non-LTE corrections are —0.3 dex 
or less severe. In the RR Lyrae case, the abundance correc- 
tions at solar metallicity are very large and ^ — 1 dex. 

At metallicities —2 < [Fe/H] < 0, the line source func- 
tion roughly follows a two-level approximation and thus the 
non-LTE corrections will depend on line strength. This ex- 
plains their increase in atmospheric models having higher 
Toft and [Fe/H], and lower log g (see Fig.|S]). 

In a more general sense (see Fig.|Sl where we show 
our resulting non-LTE corrections for a range of atmo- 
spheric parameters across the grid explored), the abun- 
dance corrections increase in size for higher effective tem- 
perature atmospheric models due to increased photon losses 
and/or level overpopulation via collisions. Corrections are 
small for TefT< 5500 K, except for solar-metallicity giants, 
where they can reach —0.6 dex. The trend with log g 
and [Fe/H] is more complex. For higher metallicity mod- 
els, the non-LTE line source function drop due to photon 
losses has a two-level nature, therefore naturally increasing 
for stronger lines. Thus, corrections tend to become larger 
with increasing temperature and metallicity/oxygen abun- 
dance and with decreasing log g in this higher metallicity 
range. Below [Fe/H]~ —2 however, the influence of log g 
and [Fe/H] is the opposite. The non-LTE corrections tend 
to increase strongly towards very low metallicity and in- 
creasing gravity, suggesting that line formation is controlled 
by a distinct non-LTE mechanism compared to the higher 
metallicity regime. 

Even though the 777 nm triplet lines get much weaker 
towards low metallicity, departures from LTE tend to in- 
crease, becoming very large for typical metal-poor turn- 
off stars. An increase in gravity at low metallicity tends 
to make the abundance corrections larger, by increasing 
the coUisional rate in the important intersystem coupling 
channel. As mentioned in Sect. I3.1?2l the non-LTE effect is 
controlled by Si i photoionization in the continuum. Thus, 
when the other parameters (Toff, log g, ^ and log eo) in the 
atmospheric model are kept constant, decreasing metallic- 
ity alone results in increasing non-LTE effects because of 
smaller continuum opacity and generally less efficient colli- 
sions (smaller contribution to free electron pool) . For exam- 
ple, even when adopting efficient H collisions (Sh = 1), the 
abundance corrections jump from ^ —0.5 to ~ —0.85 dex 
in a model with Tefr= 6500 K, log g= 4 and log eo = 6.5, 
when moving from [Fe/H]= —3 to —3.5. Finally, oxygen 
abundance, while controlling the strength of the lines, does 
not have a crucial influence on the abundance corrections, 
with differences in A log eo remaining within 0.1 dex for 
a change of 0.5 dex in log eo. 

It is appropriate to note that some authors have suggested 
to calibrate H collisions on such stars (Clementini et al. 1995; 
Gratton et al. 1999). Even though RR Lyrae may be prone to 
very large non-LTE effects, if adopting this strategy, it appears 
from our results that it would be best to calibrate on metal-poor 
turn-off stars, where the sensitivity to the choice of H collisions 
is larger 
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Fig. 8. Non-LTE abundance corrections versus metallicity, for the strongest O i IR 777 nm triplet line. Top: Tcff— 5500 K; 
log g= 2 (left panel) and log g= 4 (right panel) respectively. Bottom: same, but for Tcff = 6500 K. Results shown in all four 
panels are computed neglecting collisions with H atoms. The non-LTE corrections tend to get large at high temperature, 
both in dwarfs and giants. They are significant also in cooler solar-metallicity giants 



In summary, the non-LTE effects are particularly sig- 
nificant at the highest temperatures, in particular for low- 
gravity, solar-metallicity, and for high-gravity, very low 
metallicity models. In both cases, when H collisions are 
neglected, the non-LTE abundance corrections reach up to 
I A log eo| ~ 1 dex. The inclusion of collisions with hydro- 
gen atoms between all energy levels generally reduces the 
effects substantially (several tenths of a dex) at low metal- 
licity, while at higher metallicity it has a smaller impact, 
due to radiative losses in the line itself being mostly im- 
portant. The effects amount to —0.3 dex or less in the Sun, 
but are —0.4 dex or larger for hotter, solar-metallicity stars 
like Procyon. 

In addition, we have tested the influence of using 
ATLAS9 model atmospheres from Castelh & Kurucz (2004) 
without "convective overshooting" , instead of MARCS 
models. The temperature stratification in the MARCS and 
ATLAS models is generally thought to be quite similar 
in this parameter space (Gustafsson et al. 2008), so that 
one would expect the resulting non-LTE corrections to 
be not too different in magnitude. Any uncertainty from 
such supposedly small differences should be much less than 
that due to, say, choice of H collision efficiency or use 
of the ID approximation. We carried out a number of 
test calculations which show that, when using Castelli & 
Kurucz atmospheric models describing the Sun and stars 



of intermediately-low metallicity, the resulting non-LTE 
abundance corrections remain in fact very similar (within 
a few hundredths of a dex) to those obtained using corre- 
sponding MARCS models. The results for the most metal- 
poor cases reveal that for [Fe/H]< —2.5, the non-LTE 
abundance corrections obtained at higher-gravity become 
more significant with decreasing metallicity also for ATLAS 
model atmospheres. Table [2] lists the results obtained with 
the two sets of models. Compared to the MARCS case, the 
ATLAS non-LTE corrections, although still very signifi- 
cant, tend to be less dramatic. The largest discrepancy be- 
tween results using metal-poor MARCS and ATLAS mod- 
els is found at low gravity. For both sets of atmospheric 
models, an LTE description of the O i IR triplet is partic- 
ularly unrealistic for the lowest-metallicity ([Fe/II]= —3.5) 
dwarf case at Tcff = 6500 K. 

Knowledge of how the two sets of models compare al- 
lows us to include in the next section a discussion of our 
results in relation to other studies which have employed 
ATLAS model atmospheres. Since it is almost impossible 
to disentagle and remove all the different systematics, of- 
ten due to lack of sufficient information in previous works, 
we have chosen to compare at face value with other non- 
LTE studies in the literature, and then to apply our results 
directly to relevant existing observations. 
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Table 2. Comparison of results obtained with represen- 
tative cases from two different model atmosphere sets 
(ATLAS, and MARCS), assuming a solar oxygen abun- 
dance of log eo = 8.66 (Asphmd ct al. 2004). Wc hst the 
non-LTE corrections (obtained using our 54-level atomic 
model, respectively excluding and including the contribu- 
tion of H collisions, i.e. Sh = 0/1) for the central line in 
the Oi IR triplet. The adopted a-element abundance for 
all atmospheric models below solar metallicity is enhanced 
by +0.4 dcx 



Teff(K)\logfl 


5500\4.5 


6500\2.0 


6500\4.0 






A log eo 


A log eo 


A log eo 






(ATLAS) 


(ATLAS) 


(ATLAS) 


[Fe/Hj = 


0.0 


-0.20/-0.11 


-0.83/-0.86 


-0.53/-0.45 


[Fe/H] = 


-1.0 


-0.14/-0.05 


-0.64/-0.67 


-0.36/-0.31 


[Fe/H] = 


-2.0 


-0.il/-0.03 


-0.32/-0.35 


-0.20/-0.15 


[Fe/H] = 


-2.5 


-0.12/-0.03 


-0.28/-0.26 


-0.27/-0.16 


[Fe/H] = 


-3.5 


-0.43/-0.10 


-0.26/-0.18 


-0.87/-0.54 






A log eo 


A log eo 


A log eo 






(MARCS) 


(MARCS) 


(MARCS) 


[Fc/H] = 


0.0 


-0.18/-0.10 


-0.88/-0.92 


-0.50/-0.44 


[Fc/H] = 


-1.0 


-0.12/-0.05 


-0.63/-0.66 


-0.34/-0.29 


[Fc/H] = 


-2.0 


-0.il/-0.03 


-0.30/-0.31 


-0.19/-0.14 


[Fe/H] = 


-2.5 


-0.16/-0.04 


-0.42/-0.33 


-0.56/-0.34 


[Fe/H] = 


-3.5 


-0.52/-0.14 


-0.88/-0.53 


-1.21/-0.85 



A more detailed study comparing oxygen abundance 
corrections obtained with different sets of atmospheric 

models, should however be carried out, in particular at 
very low metallicity, where our tests reveal that when us- 
ing Castelli & Kurucz models, the non-LTE effects are less 
severe compared to results obtained using MARCS models. 
Within the parameter space explored for the two sets of 
models, the largest differences (several tenths of a dcx when 
neglecting H collisions) between the O i 777 nm triplet non- 
LTE corrections appear at the lowest metallicities, both for 
giant and dwarf models. 

The smaller departure coefficients in the very low- 
metallicity Castelli & Kurucz models is caused by reduced 
intersystem coupling (lower rates) via electron collisions, 
and thus less severe atomic level overpopulation. This is 
likely due to the effect of existing differences in electron 
density/temperature stratification in deep atmospheric lay- 
ers (for example, up to around 400 K at [Fe/H]= —2.5), 
between the MARCS and ATLAS sets of very metal-poor 
models. In this context, the failure of the Kurucz mod- 
els to provide matching oxygen abundances from the near- 
IR triplet and the forbidden line, has previously been at- 
tributed to inconsistencies in electron density and tempera- 
ture gradients of those models near the continuum-forming 
layers (Israelian et al. 2004). Compared to the MARCS 
models, we find that the effect of increasingly severe non 
LTE corrections with decreasing metallicity starts to ap- 
pear at lower metallicity in the Castelli & Kurucz mod- 
els (for example, reaching Alogeo ^ —0.7/ — 0.9 dex 
with/without H collisions, for a Tcs= 6500 K, log g= 4 and 
[Fe/H]— 4 model). We can only speculate that the two sets 
of models behave differently due to differences in opacities, 
and/or equation of state, and/or convection treatment. The 
issue of which set of ID models is more suitable or realis- 
tic is still open (see discussion in Gustafsson et al. 2008) 



and of course, current further improvements in the mod- 
elling (e.g. use of three-dimensional model atmospheres) 
may prove crucial. 

We here warn the reader that there will be system- 
atic errors when applying, as is often done, corrections 
computed for one set of model atmospheres, to LTE re- 
sults obtained with a different set of models, for example 
due to differences between continuous opacities and ion- 
ization balances used in the construction of the model at- 
mospheres and those employed for the non-LTE calcula- 
tions. Therefore, the abundance corrections we derived in 
this paper should, strictly speaking, be applied only to re- 
sults obtained with MARCS model atmospheres. However, 
given the above discussion, we suggest that they can also 
be safely applied to LTE abundances from Kurucz models 
down to intermediately-low metallicity. 



4. Comparison with other theoretical non-LTE 
studies for oxygen 

We now discuss our results in relation to those of previous 
work present in the literature. 

4.1. Kiselman 1991, 1993, 2001; and related works 

Kiselman (1991, 1993, 2001) has studied and reviewed the 
non-LTE effects on oxygen lines. A comparison with those 
works, and with CJ93, Nissen et al. (2002), Akerman et 
al. (2004), and Garcia Perez et al. (2006) is of significance. 
Our results are also relevant to the publications by Israelian 
et al. (2004), Shchukina, Trujillo Bueno & Asplund (2005) 
and Melendez et al. (2006) , where the various authors used 
atomic data based on CJ93. In fact, all the works mentioned 
in this subsection are at least in part based on the (14- and 
45-level) atomic models in that study. 

The CJ93 atomic models include estimates by van 
Rcgemortcr (1962) and the impact approximation of Seaton 
(1962) for collisions involving radiatively allowed bound- 
bound transitions, and data available from various au- 
thors for radiatively forbidden transitions connecting to the 
ground state and to the continuum. Collisional coupling 
between the singlet, triplet and quintet systems is largely 
absent in those models due to lack of data for several tran- 
sitions, in particular, for intersystem coupling via the ra- 
diatively forbidden 3s ^S° - 3s ^S° transition. Even though, 
as in CJ93, we too use the code MULTI, it is a more recent 
implementation of it (version 2.3) and we have made sev- 
eral other improvements in this work. Our atomic model 
includes a larger number of energy levels. More crucially 
with respect to the 777 nm triplet non-LTE corrections, it 
features updated atomic data, and takes into account back- 
ground opacities from lines of other elements in the region 
of the O I UV lines. 

Kiselman (1991) studied non-LTE effects on oxygen 
abundances in solar- type stars, finding that abundance cor- 
rections are slightly decreasing from solar to [Fe/H]= — 1 
but increasingly large below [Fe/H]'--^ —2, reaching as much 

as ~ -0.9 dex below [Fe/H] 2.5 (see Fig. 3 in that 

work). Fig. 4 in Kiselman (1991) shows that for a very 
metal-poor stellar model ([Fc/H]^ —3.5 dcx) the correc- 
tions found were even more severe, with a magnitude of 
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I A log eo| ^ 1-4 dex. These results, in particular of non- 
LTE effects which gently decrease down to intermediately 
low metallicity and then steeply increase at very low metal- 
licity, are qualitatively similar to those we obtain when ne- 
glecting H collisions, as Kiselman did. This is somehow for- 
tuitous due to the opposite effect on the abundance cor- 
rections at low [Fe/H] of, on one side, our newly adopted 
electron collision data, and on the other, of both our inclu- 
sion of UV background opacityand of the increased com- 
pleteness of our atomic modelEj These influences will not 
matter much for the Sun, where we too find that processes 
in the line itself are important, due to the line source func- 
tion drop having a close to two-level character. 

Our initial tests with the atomic model by K93 made 
us discover that intersystem electron collisions are crucial 
in the metal-poor regime. The original model in that work 
gives underestimated abundance corrections due to absence 
of intersystem coupling. K93 found that the introduction 
of intersystem collisional coupling changed the equivalent 
width of the 777 nm triplet by a few percent and should 
therefore be of minor importance. We confirm this if test- 
ing on atmospheric models with parameters as in that study 
(for Sun, Procyon and the metal-poor dwarf HD 140283). 
However, our calculations on atmospheric models with both 
higher temperature and lower metallicity (typical parame- 
ters of very metal-poor turn-off stars) reveal the existence 
in that range of a larger influence from intersystem cou- 
pling on the non-LTE formation of the O i IR triplet. This 
is due to the fact that in these model atmospheres there is 
an increased overpopulation of level 3s "^S", which can eas- 
ily transmit to the (common) lower level of the triplet lines 
due to its similar excitation energy. 

For the Sun, our investigation obtains non-LTE abun- 
dance corrections of ~ —0.25 dex for Sh = 0, in reasonable 
agreement with the original results of K93. We derive a cor- 
rection of ~ —0.15 dex when including inelastic H collisions 
according to the Drawin recipe. Compared to K93 (see also 
Kiselman 2001), we find that the lower level of the 777 nm 
triplet gets more overpopulated and that its upper levels 
get less underpopulated. This effect was already noted via 
test calculations in Kiselman's work. It is due to atom com- 
pleteness (i.e. to the combined infiuence of high-excitation 
energy levels) caused by the significant effect of recombi- 
nation cascades due to photon losses in high-lying infrared 
transitions. Coupling with the large population reservoir in 
the first ionization stage exists via photon suction through 
bound-bound transitions. The inclusion of high levels en- 
sures that the consequent efficient downward ffow of elec- 
trons can be well described. The balance between the con- 
tribution of source function drop and line opacity will be 
better described by our newly employed atomic model. We 
find increased contribution of line opacity. Using our atomic 
model at very low metallicity ([Fe/H]< —3), we find, in 
particular at higher effective temperature, that coupling 
between high-excitation levels and the continuum works in 
the opposite sense than in the Sun, causing an upward flow. 

Nissen et al. (2002) derived 777 nm triplet non-LTE cor- 
rections for the subset among their sample of metal-poor 
main sequence and subgiant stars in which those lines had 
been included in the observational setup. Using a previous 



Our tests show collisional coupling of high-lying levels to 
the continuum to be a significant line weakening channel at low 
metallicity 



version of our same code and the model atom in K93, and 
neglecting H collisions, they found moderate non-LTE cor- 
rections. The star with lowest metallicity (and also high- 
est temperature) in their sample is LP815-43, for which 
they derive the largest non-LTE correction of A log eo = 
—0.25 dex. When neglecting H collisions as they do, we ob- 
tain A log eo — —0.55 dex. Note however that this is likely 
an overestimate and would bring the non-LTE determina- 
tion of [0/Fe] to very low values. If including H collisions, 
the abundance correction becomes A log eo ~ —0.4 dex. 
It is interesting to note that, while their non-LTE cor- 
rections show a slight decrease with metallicity down to 
[Fe/Hj^ —2.5, there is a hint of significant increase in the 
effects at or below the metallicity of LP815-43, as stated 
at the end of their Section 4.3. This is the only star falling 
in the parameter range where we observe the very large 
increase in non-LTE corrections. Qualitatively, the trend 
with metallicity in the non-LTE corrections of Nissen and 
collaborators, reproduces our finding at overlapping metal- 
licities. Since we find collisional coupling not to play a major 
role until very low metallicity is reached, this explains why 
only the result for LPS 15-43 is significantly discrepant. The 
differences in the atomic models (namely, electron and H 
collisions, atom completeness, oscillator strengths) used in 
the two studies will indeed mostly matter in that regime. 
The results of Nissen et al. were used by Akerman et al. 
(2004) by interpolating the non-LTE corrections for rel- 
evant parameters, while Garcia Perez et al. (2006) per- 
formed tailored calculations in the same fashion as Nissen 
et al. The above conclusion - namely, underestimated non- 
LTE effects at very low metallicity - applies in those cases 
too. For all stars in the study by Garcia Perez et al., we 
obtain (neglecting H collisions as they do) more severe 
non-LTE abundance corrections, however remaining within 
^ 0.05 dex of their estimates. This would help bring their 
derived [0/Fe] values from the IR triplet in better agree- 
ment with those they determined from [Oi], however still 
without reaching full consistency between the abundance 
estimates from the two indicators. 

Shchukina, Trujillo Bueno & Asplund (2005) neglected 
H collisions in their study of the impact of non-LTE and 
granulation effects on Fe and O in the metal-poor halo 
subgiant HD 140283. Their resulting ID non-LTE cor- 
rections are smaller (—0.18 dex) than the correction they 
estimated for the Sun (which is similar to our value for 
it). Such trend agrees with our result using an appropri- 
ate model atmosphere (rcff= 5690 K, log g= 3.67 and 
[Fe/H]— 2.4), for which we too find a smaller abundance 
correction (—0.12 dex) than our estimate for a solar model. 
Note that HD 140283 lies just outside the (higher effec- 
tive temperature and lower metallicity) regime where we 
encounter large non-LTE effects. Moreover, Shchukina et 
al.'s reanalysis of its stellar parameters, aimed at solving 
the inconsistencies in the abundance determination, arrive 
at the conclusion that this star has significantly lower tem- 
perature and higher metallicity than previously thought. 



Finally, regarding comparison with the work of 
Melendez et al. (2006), they computed non-LTE correc- 
tions with the NATAJA code (see Shchukina & Trujillo 
Bueno 2001), finding a roughly constant oxygen-to-iron 
overabundance (mean [0/Fe]~ -1-0.5) from published data 
on unevolved metal-poor halo stars. Note that their Fig. 4 
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Fig. 9. Non-LTE abundance corrections versus metallicity, 
for the O I IR 777 nm triplet. Open diamonds indicate the 
results of Melendez et al. (2006), while the other points 
mark our corrections for Sh = and Sh = 1 (open and 
filled triangles respectively) 

shows that the lowest-metallicity stars tend to have smaller 
[0/Fe] values than what found around [Fe/H]^ —2. This 
may indicate (see Nissen et al. 2007) too high tempera- 
tures at low [Fe/H], a possibility either related to uncer- 
tainties in the reddening or to the not-yet settled Toff- 
scale itself. In any case, collisions with H were not included 
in their work (Melendez, private communication). A com- 
parison of the non-LTE corrections shows that theirs are 
significantly smaller at low [Fe/H]. Inspecting their Fig. 
5, where they compare with the non-LTE corrections of 
Akerman et al (2004), the star with the largest non-LTE ef- 
fect (A log eo ~ —0.35 dex) is G64-12, which has the lowest 
metallicity in their literature sample. For a model with sim- 
ilar parameters as those used by them, we obtain a larger 
abundance correction (~ —0.5 dex, or more if neglecting 
H collisions as they did). As seen in Fig.[9l while also the 
non-LTE corrections of Melendez and collaborators tend to 
become more significant at [Fe/H]< —2, we find that the 
trend with decreasing metallicity is much more pronounced, 
with more significant non-LTE abundance corrections than 
found by those authors, thus likely destroying their claimed 
quasi-fiat [0/Fe] trend with metallicity, even when adopt- 
ing the results we obtained in the case of efficient H col- 
lisions. The differences are likely related to the use of dif- 
ferent codes, model atmospheres, opacities, and to our in- 
clusion in the oxygen model atom of new electron-impact 
data. 

4.2. Gratton et al 1999 

Using an earlier version of the same code employed here 
and a similar model atom to the one in K93, but Kurucz 
(1992) model atmospheres, Gratton et al. (1999) studied 
departures from LTE in F-K stars for a broad range of 
gravities and metallicities, and including collision with elec- 
trons and hydrogen. For the latter process, they adopted 
very large efficiency (scaling factor Sh ~ 3.2 in our nota- 
tion), based on an empirical calibration on RR Lyrae stars, 
i.e. on the requirement that the 777 nm and 616 nm triplets 
should provide the same abundances. Forcing such agree- 
ment may however mean that other assumptions arc folded 



into (and remain hidden in) the scaling parameter for H 
collisions. Based on previous evidence from solar investiga- 
tions (AUende Prieto, Asplund, & Fabiani Bendicho 2004), 
there is some indication that the use of Drawin's recipe with 
no scaling factor (i.e., Sh = 1) may be reasonable for the O i 
IR triplet. Since collisions with neutral H atoms are gener- 
ally considered the most important thermalizing process in 
late- type stars, Gratton et al.'s higher choice for Sh will give 
smaller non-LTE corrections than in our case. Indeed, those 
authors found relatively mild non-LTE effects on the Oi 
777 nm triplet, in particular —0.1 dex in the Sun. Their 
study reaches down to [Fe/H]= —3, where for hot, high- 
gravity atmospheric models they find the largest correc- 
tions (reaching A log eo ^ —0.5 dex) for the Oi IR triplet. 
Choosing atmospheric parameters in common between the 
two studies {Tcs= 6000 K, log g= 4.5, [Fe/H]= -3, see 
their Fig. 12), and very efficient H collisions as they did, 
we too obtain small non-LTE abundance corrections (only 
^ —0.07 dex in our case, i.e. even a few hundredths of a 
dex less severe than in Gratton et al.). However, if adopting 
Sh < 1, our non-LTE corrections will actually be larger 
(A log eo < —0.25 dex). Differences with our study are 
therefore mostly due to our use of H collisions more in 
line with the majority of presently available evidence, and 
only to a lesser extent to employing different model atmo- 
spheres. It is interesting to note (their Fig. 12) that their 
non-LTE corrections for the hottest model used by them 
(Teff= 7000 K) show a significant increase when going from 
[Fe/H]= —2 to —3. The results of Gratton and co-workers 
were for example included in Garretta, Gratton & Sneden 
(2000) and in Primas et al. (2001), deriving a "quasi-flat" 
oxygen-to-iron overabundance with mean [0/Fe]'^ -1-0.5 at 
low metallicity. 

4.3. Takeda 1992, 1994, 2003; Takeda & Honda 2005 

Takeda published a series of non-LTE studies of oxygen. 
In his investigations, rates due to collisions with H atoms 
were in general computed using the Drawin recipe without 
corrections. 

Takeda (1994) found [Alog eo| < 0.1 dex for the 
Sun but ~ —0.4 dex for a hotter model representative of 
Procyon (for which, when including H collisions, we derive 
a correction of— 0.4/ — 0.5 dex depending on the O i triplet 
line considered). That author also attributed to non-LTE 
effects, at least qualitatively, the discrepancy between oxy- 
gen abundances from the triplet and forbidden lines gener- 
ally found in studies of low- [Fe/H] halo stars. In Takeda's 
non-LTE analysis, the most oxygen-deficient (log 0= 6.86) 
calculation starts to show the signs of UV pumping in the 
resonance lines (see their Fig. 9 and discussion in the cor- 
responding section). 

Takeda (2003) gives an analytical formula, indepen- 
dent of metallicity, to approximate the non-LTE abun- 
dance corrections resulting from his calculations, arguing 
that it will provide a fairly good approximation for lines 
with Wa < 100 mA. For the Sun (see also Takeda & Honda 
2005) and for solar-type dwarfs at moderately low [Fc/H], 
our calculations including H collisions give similar results 
to those of Takeda. The mentioned formula by Takeda does 
not explicitly contain [Fe/H]. However, given its dependence 
on line strength, it will give generally less severe abundance 
corrections at low metallicity due to weaker lines, so that 
our non-LTE corrections become larger than Takeda's at 
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low metallicity. Note that his Fig. 4 shows that estimated 
corrections are low and roughly constant —0.1 dex) for 
the stars considered in his comparison with Nissen et al. 
(2002). As discussed in 14. 11 we obtain roughly similar non- 
LTE corrections to Nissen et al. (2002), or even more se- 
vere effects for the lowest-metallicity star in their sample. 
Note however that when adopting efficient H collisions as 
in the study by Takeda, the results will be closer. The ab- 
sence of a significant trend of the non-LTE corrections with 
decreasing metallicity found by Takeda (2003) is at vari- 
ance with the conclusions of e.g. Carretta et al. (2000). 
Takeda argues that, even though in shallow atmospheric 
layers there is an appreciable increase of the non-LTE line 
opacity at low metallicity, this does not affect the deeply- 
forming 777 nm O i triplet, so that fine details of the atomic 
model will not matter, due to a close-to-two-level line source 
function. Note that Takeda does in fact find a metallicity 
dependence when neglecting H collisions. We find that at 
low metallicity, the two-level nature of the non-LTE effect 
breaks down. In particular, for Tcs= 6500 K, log g= 4, us- 
ing Takeda's formula with our estimated Wnon-LTE would 
give a constant A log eo ^ —0.1 dex at very low metallic- 
ity, compared to our results of ^ —0.5 dex and ~ — 0.9 dex 
(respectively for [Fe/H]= —3 and —3.5) when adopting H 
collisions. In this respect, our abundance corrections would 
be more helpful in solving the oxygen discrepancy, bringing 
the abundances derived from the 777 nm triplet in bet- 
ter agreement with the usually lower value from [O i]-based 
LTE estimates. 

Thus, we do not confirm the small dependence of Oi 
non-LTE corrections from [Fe/H] found by Takeda. That 
author (Takeda, private communication) included intersys- 
tem collisional couphng using Auer & Mihalas (1973) and, 
in Takeda (1992) and successive works, treatment of scat- 
tering in the continuum source function. Note that Takeda 
accounted for background line opacity in the UV only from 
H lines. However, as seen from our tests, metal lines tend to 
make the non-LTE effect smaller. Therefore, this actually 
makes the discrepancy with our results more significant, 
and likely explainable due to the summed effect of different 
atomic data, and use of different model atmospheres. The 
results in Takeda (2003) were applied by means of interpo- 
lation/extrapolation in Takeda & Honda (2005), giving a 
linear increase of [0/Fe] with [Fe/H]. 

4.4. Ramirez, Allende Prieto & Lambert 2007 

In their study of a large sample of nearby stars, Ramirez, 
Allende Prieto & Lambert (2007) have derived non-LTE 
corrections (neglecting the role of inelastic H collisions) for 
a grid of model atmospheres, mainly covering the metal- 
licity regime of the Galactic thin and thick disks. Their 
adopted O i model atom is explained in detail in Allende 
Prieto et al. (2003). It consists of 54 levels plus continuum 
and 242 transitions, with the radiative data coming from 
the Opacity Project and NIST, and electron collision rates 
estimated via approximate formulae. For excitation by elec- 
trons in forbidden transitions, they computed the collision 
rate coefficient from Eissner & Seaton (1974), adopting an 
arbitrary choice of effective collision strength. As shown in 
the detailed comparison by Barklem (2007a) , the data here 
employed by us for radiatively forbidden transitions show 
poor agreement with those estimated by Allende Prieto et 
al. In particular, for the radiatively forbidden transition 



which we find important in coupling triplet and quintet sys- 
tems, our adopted rate is two orders of magnitude larger. 

The non-LTE corrections of Ramirez et al. reach down 
to - -0.4 dex at [Fe/H]-- -1.5 (their Fig. 7). For stars 
with determined metallicities, as from their Table 6, the 
largest corrections are ~ —0.3 dex (e.g., HD 112887 at 
Teff= 6319 K, log g= 3.88, [Fe/H]== -0.38). Their abun- 
dance corrections show a relatively mild dependence on 
metallicity for the range studied. Since our investigation 
extends to much lower metallicity than theirs, we only com- 
pared the results for a few representative stars with atmo- 
spheric parameters covered by both works, using the rou- 
tine (Ramirez, private communication) mentioned in their 
paper. The star with lowest determined metallicity in Table 
6 of Ramirez et al. (2007) has [Fe/H]= -1.46. For this star 
(HD 132475), their routine gives a non-LTE correction of 
~ —0.28 dex with their adopted parameters of rcff= 5613 K 
and log g= 3.81. Our non-LTE correction for similar pa- 
rameters is smaller, ^ —0.17 dex. 

For a standard star like Procyon (rcff= 6530 K, log g= 
3.96, [Fe/H]= 0), for which their formula would give ^ 
—0.28 dex, we derive a much larger estimate of ^ —0.55 dex, 
when adopting Sh — 0. Note however that, as discussed, a 
choice of Sh ^ 1 niay be more appropriate, which would 
lead to a value of A log eo ~ —0.45 dex in our case. 

For the Sun, their non-LTE correction is roughly ^ 
—0.15 dex, which is '^0.1 dex less severe than in our case. 
Even adopting their slightly higher (LTE) solar oxygen 
abundance, our results for the Sun would not change notice- 
ably. A choice of Sh ^ 1 for H collisions in our case would 
bring, however fortuitously, our estimate close to that of 
Ramirez et al. 



5. Implications 

Our results have important consequences in relation to 
the debate regarding the [0/Fe] trend at low metallic- 
ity, and on the understanding of the chemical evolution of 
this element in the early Galaxy. The large difference of- 
ten found in the literature between the oxygen content of 
old stars as derived via the O i IR triplet and the other 
oxygen abundance indicators (in particular, [Oi]) is likely 
to be explained in terms of various uncertainties acting to 
produce disagreeing estimates (Asplund 2005; Melendez et 
al. 2006). 

In our re-assessment of non-LTE effects on the O i IR 
triplet, we have found that at low metallicity they are more 
important than what is often assumed in the literature. 
In Fig. 1101 we show the effect of applying our abundance 
corrections to literature data by Israelian et al. (2001), 
Nissen et al. (2002) and Melendez et al. (2006). Since un- 
certainties in H collision efficiency affect the modelling of 
the 777 nm O i triplet (and corresponding non-LTE cor- 
rections) very significantly, we have plotted results derived 
applying Sh = and Sh = 1, respectively. The limited ob- 
servational evidence for the Sun suggests for oxygen rather 
efficient thermalization via impacts with H atoms, akin 
to choosing Drawin's formula without any scaling factor 
(Sh = f )• The severe non-LTE corrections we find seem to 
largely erase the linear increasing trend towards low metal- 
licity sometimes found in LTE. Our estimates with a choice 
of Sh = 1 suggest a rather fiat trend below [Fe/H]= —1, al- 
though with significant residual scatter. Only adopting an 
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Fig. 10. Revised [0/Fe] versus [Fe/H], from literature data. The different symbols are: filled circles, for [OI]-based, 
LTE-obeying values from Nissen et al. (2002); filled triangles, for the ID estimates (from that same study, but non-LTE 
corrected according to our results) based on the Oi 777 nm triplet; filled diamonds for Oi-based values for the stars 
in Israelian et al. (2001) with fully determined stellar parameters and oxygen triplet abundances in their study, taking 
into account our non-LTE corrections on the oxygen triplet; and empty squares, for our non-LTE corrected estimates 
using the Melendez et al. (2006) data set (after adjusting Tcg-scale, see Nissen et al. 2007). The non-LTE abundance 
corrections we applied were computed with a choice of Sh = and Sh = 1 (left and right panel respectively) for hydrogen 
collisions. The dotted hue indicates the mean trend suggested by Israehan et al. (2001), i.e. [0/Fc]= — 0.33[Fc/H] 



unrealistically high amount of thermalization would make 
the non-LTE effects negligible. When instead completely 
ignoring H collisions, (as discussed however, this case may 
be unrealistic), the [0/Fe] values in Fig. [TO] would become 
even lower (the trend turning into a decrease with decreas- 
ing metallicity) . If such a model was borne out, it would 
possibly require new investigations on the efficiency of the 
^^C(a, 7)^^0 nuclear reaction and on published nucleosyn- 
thetic yields at low metallicity. 

We note that our non-LTE corrections do not destroy 
the good agreement between [O i] and O i abundances 
found by Nissen and collaborators in five of their sample 
stars. Since those stars are between —2.42 < [Fe/H] < 
— 1.04 they do not suffer the very large corrections (nor 
the large sensitivity to H collisions) we found at very low 
metallicity. In fact, our resulting non-LTE corrections, and 
thus estimated abundances from the 777 nm triplet, agree 
with those of Nissen et al. (2002) within 0.1 dex for those 
stars. 

Our results may imply that the oxygen-to-iron over- 
abundance does not depend significantly on metallicity 
down to [Fe/H]= —3, with an essentially flat trend (this 
may include a slight slope in either direction) at very 
low metallicity. This is in line with the behaviour recently 
found via accurate abundance estimates (taking into ac- 
count non-LTE effects) for other typical a-capture elements 
(e.g. Nissen et al. 2007). In Fabbian et al. 2008b, we further 
investigate this issue, together with a study of the [C/Fe] 
and [C/0] ratios at low metallicity, via new observational 
data. 

We have shown that O i LTE abundances derived from 
the 777 nm triplet are likely large overestimates, espe- 
cially at low [Fe/H], but there is still significant uncer- 
tainty related to the choice of H collision efficiency. Some 
authors have argued on astrophysical grounds that assum- 
ing the Drawin estimate might be reasonable, at least for 
the Sun. Solar center-to-limb variation in fact not only pro- 
vides evidence that the O i 777 nm lines are formed in 



non-LTE (Altrock 1968; Shchukina 1987; Kiselman 1991), 
but also suggests a possible significant role of H collisions 
(e.g. AUende Prieto, Asplund, & Fabiani Bendicho 2004). 
Relative to electron collisions, these should be even more 
important at low metallicity. The model with Sh = 1 seems 
to predict lines in better agreement with astrophysical ex- 
pectations from other oxygen spectral features. However, 
this is only an indication, since it could well result from er- 
rors due to the Drawin estimates and to the use of ID atmo- 
spheric models coincidentally acting in opposite directions. 
Here, we too find that H collisions may need to be taken 
into account for this element, by reasoning on the otherwise 
unlikely large impact of non-LTE corrections on [O /Fe] at 
very low metallicity, when they are neglected. At the same 
time, we warn that in fact (however unavoidable in the ab- 
sence of available detailed quantum mechanical computa- 
tion for relevant 0-l-H collisions, which are urgently needed 
together with constraints via solar observations) adopting 
the general Drawin recipe - possibly with some scaling fac- 
tor - for all transitions is a rough approximation, since dif- 
ferent transitions likely have different sensitivity to such 
collisions. 

Realistic cross-sections for continuum absorption are 
also needed. We have found that Si abundance plays an 
important role (via bound-free absorption) in the forma- 
tion of the IR oxygen triplet at very low metallicity. This 
element is commonly thought to follow the normal be- 
haviour of other a-elements in being overabundant com- 
pared to Fe at low metallicity. Cayrel et al. (2004) recently 
obtained an average [Si/Fe]~ -1-0.4 (with small scatter), 
which we have adopted in the calculations. However, note 
that very recently Nissen & Schuster (2008) have found evi- 
dence that halo stars may fall into two groups, with distinct 
a-enhancement. If this significant (~ 0.20 dex) separation 
between "low-a" and "high-a" halo stars is preserved down 
to very low [Fe/H] for silicon, it may well be an important 
factor in the discussion of the oxygen problem, because a 
lower content of Si would imply larger non-LTE effects on 
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the Ol triplet, for an otherwise fixed set of atmospheric 
parameters. Any residual star-to-star differences in the Si 
content of very metal-poor halo stars may thus play some 
role in the scatter observed in [0/Fe] versus [Fc/H]. It is 
therefore important, in future studies on the oxygen prob- 
lem, to pin down the non-LTE corrections by also having 
very accurate determinations of the stellar silicon content. 

We cannot yet completely rule out the existence of a 
continuous steeply increasing trend of [0/Fe] with decreas- 
ing [Fc/H]. However, its existence would require extremely 
efficient H collisions (e.g. Sh > 3), in order to bring the 
abundance estimates close to the high extreme set by the 
LTE expectation. 

The oxygen-to-iron overabundance more likely reaches 
^ +0.5 dex at very low metallicity, with a plateau or slight 

metallicity-dependence. Further modelling and new high- 
quality observations will be needed to completely settle 
the issue, but models of Galactic chemical evolution con- 
strained by the [0/Fe] ratio derived from observations of 
stars at low metallicity, can now hope to achieve a more 
definite conclusion on possible scenarios required by de- 
rived stellar composition at early stages. Such low (com- 
pared to LTE) [0/Fe] values in Galactic halo stars are in 
better agreement with recent estimates (Pequignot 2008) 
for low-metallicity blue compact dwarf galaxies. 

Regarding the solar oxygen abundance, it is little af- 
fected by the new electron collisions. The determination 
of Asplund et al. (2004) employed several different indi- 
cators obtaining remarkable agreement between them, by 
performing a 3D non LTE analysis. Our ID non LTE re- 
sults for the 777 nm IR triplet lines agree within ^ 0.05 dex 
with those published by Asplund et al., when adopting their 
choice of Sh = 0. We argue that, even if H collisions were 
important for oxygen (Asplund and collaborators neglected 
them in their calculations) - in which case our estimates 
would give up to ~ —0.1 dex higher oxygen abundance 
for some lines (in particular, for the 844.67 nm line) due to 
smaller corrections than theirs - the overall change when 
averaging among all lines will remain within 0.05 dex. This 
is still significant, since the derived solar oxygen abundance 
may be around ~ 8.70, which would partially alleviate the 
current discrepancy with helioseismology results. 

Incidentally, we cannot explain the extremely large 

[O /Fe] overabundances derived in some cool giants at very 
low metallicity, and the very different abundances derived 
in those objects from permitted and forbidden oxygen lines 
respectively (Israelian et al. 2004), because in that range 
our non-LTE corrections tend to become ~ —0.1 dex or less 
severe. To gain a better insight into the various outstanding 
problems, detailed results from 3D non-LTE investigations 
will be necessary. 

In a series of papers, Schuler et al. (2004, 2005, 2006) 
have studied solar-type open cluster dwarfs, finding much 
lower mean [Oi] -based abundances than for results based 
on the permitted Oi triplet. From the latter, they found 
a puzzling and abrupt increase in LTE oxygen abundance 
derived from the high-excitation triplet for stars with Teff 
< 5800 K. We exclude that this trend could be due for 
the most part to departures from LTE, since our resulting 
abundance corrections tend to be relatively small in the 
relevant temperature range. 



6. Conclusions 

We have explored the oxygen non-LTE line formation over 
a wide parameter space. Thanks to several improvements - 
mainly, detailed treatment of UV radiation field including 
opacities by background lines, and extended model includ- 
ing recently available atomic data - we derive estimated 
corrections at low [Fe/H] to LTE abundances derived from 
the IR 777 nm triplet that are generally larger than usually 
adopted in the literature. It is clear that the availability of 
adequate atomic data is crucial to achieve a reliable non- 
LTE solution for this and other important astrophysical 
problems. Using rates derived from newly available esti- 
mates for electron-impact excitation by electrons (Barklem 
2007a), we have performed non-LTE calculations in one- 
dimensional stellar atmospheric models of late-type stars 
extending to very low metallicity, in order to understand 
the formation of the permitted neutral atomic oxygen lines. 
The trends of non-LTE corrections across the parameter 
space can be understood in terms of the different non LTE 
mechanisms in action: mainly source function drop at solar 
metallicity, and large level overpopulation at low [Fe/H]. 

For solar-type dwarfs, the non-LTE corrections we 
found show relatively small metallicity-dependence down 
to [Fe/H]~ —2. This is in line with previous findings by 
other authors (e.g., Takeda 2003). However, at the higher 
temperatures of turn-off stars at low [Fe/H], the departures 
from LTE tend to increase to large values when metallicity 
decreases, so that LTE abundances will require increasingly 
large corrections. The two- level approximation valid for the 
Sun breaks down and results no longer mostly depend on 
processes in the line itself. The large non LTE corrections 
at low [Fe/H] follow from increased coUisional intersystem 
coupling. The (over)population of the lower level of the 
777 nm triplet tends to increase and so does line opacity. 
The formation of the 777 nm lines becomes sensitive, at low 
metallicity, to absorption processes in the UV continuum. 
This makes the non-LTE corrections metallicity-dependent 
due to the influence of Si absorption. The results have an 
impact on the derivation of the "true" [0/Fe] trend at low 
metallicity, with a roughly constant plateau from literature 
data if including thermalization by collisions with hydrogen 
atoms. 

The problem of H collisions is obviously still very much 
open. Progress is required, via theoretical quantum me- 
chanical calculations or via limb darkening observations 
(AUende Prieto, Asplund & Fabiani Bendicho 2004), in 
order to assess whether H collisions can provide efficient 
thermalization. This is of high priority also for a more con- 
clusive determination of the solar oxygen content as stan- 
dard of reference in abundance determinations. Some au- 
thors (e.g., Kiselman 1991; Nissen et al. 2002; Asplund et 
al. 2004) have preferred to completely neglect H collisions, 
because of experimental (e.g. Fleck et al. 1991) and the- 
oretical (e.g. Barklem, Belyaev & Asplund 2003) evidence 
that Drawin's formula gives usually too large (by orders of 
magnitude) estimates for atmospheres of cool stars, at least 
for simple atoms like Li. We find that even with a choice of 
Sh = 1, the non-LTE corrections would still remain large 
at very low metalHcity (up to —0.85 dex at [Fe/H]= —3.5). 

Even though when adopting extreme choices of the pa- 
rameter regulating the efficiency of H collisions, the differ- 
ence in the resulting 777 nm triplet non LTE corrections is 
relatively small towards higher metallicity, it remains sig- 
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nificant and '--^ 0.1 dex for the Sun, in the sense of efficient 
collisions producing smaller non-LTE effects (which would 
help to partially alleviate the discrepancy with helioseis- 
mology, at least for oxygen). 

It looks as though a combination of 3D/non-LTE ef- 
fects, choice of atmospheric models and of temperature 
scale, and observational uncertainties both in derived stel- 
lar parameters and equivalent widths measurements have 
conspired to prevent a clear solution of the oxygen problem 
so far. It now seems like better agreement with the abun- 
dance from the forbidden [O i] lines (which should be free 
from non-LTE effects) and from OH molecular lines (once 
corrected for 3D effects, Asplund & Garcia Perez 2001) can 
be obtained when non-LTE effects are taken into account 
for the 777 nm oxygen triplet abundance, see e.g. Melendez 
et al. (2006). Our results seem to suggest that large non- 
LTE corrections at low metallicity are the norm, at least 
for turn-off halo stars. 

Given our findings concerning the importance of the 
radiatively forbidden 3s '^S" - 3s transition, it is crucial 
to check if a similar behaviour is seen for other atoms with 
similar structure, e.g. C, N, and S, or even in all atoms. 

Non-LTE corrections obtained with ATLAS models 
turned out to be very similar, for a range of atmospheric 
parameters, to those we found using MARCS models. 
However, at very low metallicity, significant discrepancies 
appear, with large differences between non-LTE corrections 
using the two sets of models. Interestingly, this metallic- 
ity range overlaps with that for which the debate between 
"flat" and "linear" [0/Fc] trend exists. However, as seen 
from Fig. 1101 our main conclusions will not change even if 
applying ATLAS non-LTE corrections, which would give, 
around [Fe/H]= —2.5, smaller corrections by ~ 0.3 dex in 
the left panel (no H collisions) and by 0.2 dex in the right 
panel (H collisions 'a la Drawin'). Anyhow, more sophisti- 
cated modelling of stellar atmospheres is likely another cru- 
cial factor with respect to the oxygen problem, in particular 
for abundances from OH lines, but also to understand how 
the large non-LTE effects we found act in a 3D atmosphere. 
Therefore, it is obviously urgent to carry out, for a large 
range of stellar parameters, full 3D non-LTE calculations, 
in order to further investigate the formation of the oxygen 
lines at low [Fe/H]. It is expected that non-LTE effects 
will be enhanced when using 3D atmospheric models, due 
to their cooler superficial temperatures, to which the high- 
excitation lines of interest are known to be sensitive. This 
crucial step forward will allow enhanced abundance anal- 
yses and will help to shed light on early Galactic chem- 
ical evolution, contributions by SN II and massive stars, 
time delay of yields from Type la SNe, and the reality of 
the [C/0] upturn at low metallicity (Akerman et al. 2004; 
Spite et al. 2005; Fabbian et al. 2008b). It is hkely that 
- together with the use of an accurate temperature scale - 
when finally taking into account, using detailed calculations 
with the best available atomic data, the interplay of 3D and 
non-LTE effects in the formation of the various available 
abundance indicators (including for the determination of 
metallicity), it will be possible to in the end fully solve the 
"oxygen problem" . 

Here, we have demonstrated that non-LTE corrections 
play an important role towards a solution. They will need 
to be taken into consideration in order to abate large sys- 
tematic errors still afflicting LTE-based estimates. To check 
whether agreement between abundances from the different 



oxygen lines can be found at low metallicity will further- 
more require using improved observational data, in par- 
ticular from a large sample of subgiants or giants where 
more oxygen indicators are available. This important test 
will help to clarify outstanding issues related to the [O /Fe] 
controversy. 
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